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We review the method of blocking of topological defects from continuum used as 

a non-perturbative tool to construct effective actions for these defects. The actions 

- - - are formulated in the continuum limit while the couplings of these actions can be 

f— ^ . derived from simple observables calculated numerically on lattices with a finite lattice 



(N 



^D ! spacing. We demonstrate the success of the method in deriving the effective actions 

for Abelian monopoles in the pure SU(2) gauge models in an Abelian gauge. In 
particular, we discuss the gluodynamics in three and four space-time dimensions at 
zero and non-zero temperatures. Besides the action the quantities of our interest 

cn . are the monopole density, the magnetic Debye mass and the monopole condensate. 

V. PACS numbers: 11.15.Ha,14.80.Hv,11.10.Wx 
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,Q ; I. INTRODUCTION 

\D 
O 
"^ ' The blocking from continuum (BFC) is a well-known tool to construct the "perfect 

actions" for lattice field theories [l|. By definition a perfect lattice action does not depend 
c^ ! on the cut-off parameter which is usually associated with the finite lattice spacing. The cut- 

off dependence provides a systematic error in the lattice observables which is of the order 
of the lattice spacing for the standard Wilson action. The various improvement schemes [2] 
^ . are used to decrease the cut-off influence on the lattice results, and the BFC method [l| is 

one of the practically useful tools used in the lattice simulations nowadays. 

Although the main idea of introducing the BFC is to reduce the systematic errors in the 
numerical simulations, a philosophically similar method can be applied to various topolog- 
5t , ical defects. As a result one can derive effective actions for the defects in the continuum 

limit using results of the lattice simulations obtained on lattices with flnite lattice spacings. 
This short review is devoted to a demonstration of success of the method applied to the 
Abelian monopoles in the lattice gluodynamics in three dimensions (where the monopoles 
are instanton-like objects)and in four dimensions (where the monopoles are particle-like 
defects) following Refs. [s], y, |5| . One should stress from the very beginning that the method 
is quite general and is not limited to the lattice monopoles only. 

First, let us remind basics of the BFC method for fleld degrees of freedom. A simplest 
way to associate, say, a continuum free fermion fleld, ^(x), with a lattice fermion fleld, \l/s, 
isQ: 

■^s= I d^xV^(x), ^s= I d^xV5(x), (1) 
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where the integration is carried out over the lattice hypercube, C^, centered in the lattice 
point s (we will come to the precise definition of Cg later). Equations (P) can then be 
inserted into the partition function as the (5-function constraint. To complete the procedure 
of blocking the continuum fields ipi^) should be integrated out leaving us with the partition 
function depending solely on the lattice fields \E's. Similar relations can also be written for 
the gauge fields etc.. We refer a reader interested in the blocking of fields to the original 
articles switching at this point to the blocking of the topological defects. 

Suppose, that we have a (gauge) model which describes topological defects, say, for 
definiteness, monopoles. In four space-time dimensions (4D) the monopole is a particle-like 
object {i.e., its trajectory is line-like), and the monopole charge is quantized and conserved 
(i.e., the monopole trajectories are closed loops). Obviously, basic requirements to the 
topological BFC procedure should be the following: (i) the procedure should provide us 
with the configuration of the lattice monopole currents for a given configuration of the 
continuum monopole currents; (ii) the lattice monopole currents must be closed; (iii) the 
lattice magnetic charge for such monopoles must be quantized. We show below that one may 
write a blocking relation similar (but, in general, not identical) to Eq. (Q), which associates 
the lattice and the continuum monopole charges and preserves their basic properties. We 
insert this relation into the partition function in a form of the (5-function constraint, integrate 
out the continuum degrees of freedom and get the lattice model which contains only the 
lattice monopole currents. Using the BFC method one can also get analytical formulae for 
various lattice observables expressed through the parameters of the continuum model. A 
comparison of the numerical data for such observables with the corresponding analytical 
expressions provides us with the parameters of the monopole action in the continuum. Note 
that the blocking of the topological defects from the continuum to the lattice is similar to 
ideas of Refs. p, |7| which discussed theoretically the blocking of the monopoles from fine 
lattices to coarse lattices. 

Below we describe how the method works for the Abelian monopoles in SU(2) gluody- 
namics. Many monopole observables have been calculated numerically [8|. Even an almost 
perfect monopole action on the lattice has been determined (in a truncated form) using an 
inverse Monte-Carlo method J9|. However, the correctness of the truncation, or, in other 
words, the correct form of the perfect lattice action is not known. The BFC method allows 
us to find couplings of the (truncated) perfect monopole action in the continuum, estimate 
the error of the truncation, and to obtain certain non-perturbative quantities. Our interest 
in the physics of the Abelian monopoles in the non-Abelian pure gauge theories is stimu- 
lated by the relation of the monopole dynamics to the one of the most important problems of 
QCD, the confinement of color. One popular approach to this problem is the so-called dual 
superconductor mechanism llOll (for a review of another interesting approach, the vacuum 
correlator method, see Ref. |ll|). The key role in the dual superconductor mechanism is 
played b y A belian monopoles which are identified with the help of the Abelian projection 
method lj|. The basic idea behind the Abelian projections is to fix partially the non- 



Abelian gauge symmetry up to an Abehan subgroup. For SU(N) gauge theories the residual 
Abelian symmetry group is compact since the original non-Abelian group is compact as 
well. The Abelian monopoles arise naturally due to the compactness of the residual gauge 
subgroup. 

The Abelian monopoles are not present in QCD from the beginning: they are not solutions 
to the classical equation of motion of this theory. However, the monopoles may be considered 
as effective degrees of freedom which are responsible for confinement of quarks. According to 



the numerical results |13| the monopoles are condensed in the low temperature (confinement) 
phase. The condensation of the monopoles leads to formation of the chromoelectric string 
due to the dual Meissner effect. As a result the fundamental sources of chromoelectric field, 
quarks, are confined by the string. The importance of the Abelian monopoles is also stressed 
by the existence |l^ of the Abelian dominance phenomena which were first observed in the 
lattice SU{2) gluodynamics: the monopoles in the so-called Maximal Abelian projection |l5l | 
make a dominant contribution to the zero temperature string tension. 

In the deconfinement phase (high temperatures) the monopoles are not condensed and 
the quarks are liberated. This does not mean, however, that monopoles do not play a 
role in non-perturbative physics. It is known that in the deconfinement phase the vacuum 
is dominated by static monopoles (which run along the "temperature" direction in the 
Euclidean theory) while monopoles running in spatial directions are suppressed. The static 
monopoles should contribute to the "spatial string tension" - a coefficient in front of the 
area term of large spatial Wilson loops. And, according to numerical simulations in the 
deconfinement phase |l6|, the monopoles make a dominant contribution to the spatial string 
tension. Thus, the monopoles may play an important role not only in the low temperatures 
but also in the high temperature phase. 

We refer a reader to Ref. [8| for a review on the Abelian projections and the dual super- 
conductor models in non-Abelian gauge theories. 

The structure of the paper is as follows. In Section|n]we describe the BFC method in the 
simplest three-dimensional (3D) case. Assuming that in the continuum the monopole action 
is of the Coulomb form we derive the lattice monopole action and the lattice density of the 
(squared) monopole charges. In Section IIIII we apply the BFC procedure to the Abelian 
monopoles in the four-dimensional SU(2) gauge theory. Assuming that the monopoles 
are described by the dual Ginzburg-Landau model one can get an analytical form for the 
quadratic monopole action on the lattice. 

Then in Section HVl we compare the analytical formulae with the corresponding numerical 
data obtained in the three dimensional SU(2) gauge model. As a result we get the density 
of the monopoles and the monopole contribution to the magnetic screening length in the 
continuum limit of this model. We show that the results obtained with the help of the BFC 
method are in agreement with the results of other (independent) calculations. In Section IVI 
we apply the 3D BFC method to the temporal components of the monopole currents in 
SU(2) gauge model in the four space-time dimensions at high temperature. This gives us a 
numerical value of the product of the Abelian magnetic screening mass and the monopole 
density in the continuum model. The self-consistency check shows that the dynamics of 
the static monopole currents can be described by the Coulomb gas model starting from the 
temperatures T > 2.5Tc. 

Finally, in Section IVII we get the value of the monopole condensate in the continuum 
using the 4D BFC method. This value is in agreement with the results obtained with the 
help of other methods. Our conclusion in presented in the last Section. 

II. BLOCKING IN THREE DIMENSIONS 

It is instructive to start the description of the BFC method from the simplest three- 
dimensional case. In three dimensions the Abelian monopoles are point-like objects char- 
acterized by a position, x, and the magnetic charge, q (measured in units of a fundamental 
magnetic charge, qm)- The simplest model possessing the monopoles is the 3D compact 



quantum electrodynamics (cQED3)in which the monopole action is given by the SD Coulomb 
gas model |l7| : 
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The Coulomb interaction in Eq.Q is represented by the inverse Laplacian D: —dfD{x) = 
5^^\x), and the latin indices a,b label different monopoles. To get analytical expressions 
below we make a standard assumption that the density of the monopoles is low. The 
monopole charges therefore are restricted by the condition \qa\ < 1 which means that the 
monopoles do not overlap. The average monopole density p is controlled by the fugacity 
parameter (, giving p = 2^ in the leading order of the dilute gas approximation [l7| 

The magnetic charges in the Coulomb gas (0) are screened: at large distances the two- 
point charge correlation function is exponentially suppressed, (p(x)p(y)) ~ exp{ — |x — y\/ 
Xd}- Here Ad is the Debye screening length [17|, 

Ad = ^. (3) 

The three dimensional Debye screening length corresponds to a magnetic screening in four 
dimensions. Below we choose the vacuum expectation value of the continuum monopole 
density p and the Debye screening length A^ as suitable parameters of the continuum model 
(instead of qm and Q. 

Next, let us consider a lattice with a finite lattice spacing b which is embedded in the 
continuum space-time. The cells of the lattice are defined as follows: 

Cs = lb[s^-^)<x^<b[si + ^), z = l,2,3|, (4) 

where Si is the lattice dimensionless coordinate and Xi corresponds to the continuum coor- 
dinate. 

The basic idea of the BFC method is to treat each lattice 3D cell as a "detector" of the 
magnetic charges of the continuum monopoles. The relation between the lattice magnetic 
charge kg and the density of the continuum monopoles p(x) is^ 

ks= fd'xpix), p{x) = Y.qa5^''\x-x^'^), (5) 

Cs 

see an illustration in Figure H The fluctuations of the monopole charges of the lattice 
cells must depend on the properties of the continuum monopoles. As a result, the lattice 
observables - such as the vacuum expectation value of the lattice monopole density - must 
carry information about dynamics of the continuum monopoles. The observables should 
depend not only on the size of the lattice cell, 6, but also on the features of the continuum 
model which describes the monopole dynamics. 



^ This relation is similar to the blocking of the continuum fields (^. In the four dimensions this similarity 
disappears. 
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FIG. 1: Blocking of the continuum monopoles to the lattice in three dimensions. The charge of the 
lattice monopole in the cube C is given by the total magnetic charge of the continuum monopoles 
inside this cube. 

It is worth stressing the difference between the continuum and the lattice monopoles: the 
continuum monopoles are fundamental point-like objects^ while the lattice monopoles are 
associated with the finite-sized lattice cells with non-vanishing total magnetic charge. 

According to definitions (0), the lattice monopole charge shares similar properties to 
the continuum monopole charge. The monopole charge ks is quantized, ks G ^, and it is 
conserved in the three-dimensional sense: 



s<=A i 



(6) 



if the continuum charge is conserved. Here A and V denote the lattice and the continuum 
volume occupied by the lattice, respectively. In other words, the total magnetic charge of the 
lattice monopole configuration is zero on a finite lattice with periodic boundary conditions. 
In next two subsections we follow Ref. ]^ presenting in the BFC approach the simplest 
quantities characterizing the lattice monopoles: the monopole action Smon{k) and the vac- 
uum expectation value of the squared magnetic charge, (fc^). 



A. Monopole action in 3D 

To derive the lattice monopole action we substitute the unity, 

1= E n<^A(A:.-/d3a;p(:5 



fc(A)e^seA 



(7) 
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^ In fact the Abelian SU(2) monopoles have a finite core |l^ of the order of 0.06 fm, which is neglected in 
our approach. 
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into Eq.(j2I). Here J2k{A)e^ — YiseAJ^kse^ ^^^ '^A stands here for the Kronecker symbol (i.e., 
lattice ^-function). We get 
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where we have introduced two additional integrations over the continuum field x ^i-nd the 
compact lattice field h to represent the inverse Laplacian in Eq.® and the Kronecker symbol 
in Eq.((Zj), respectively. The subscript A in V^^^h indicates that the integration is going over 
the lattice fields h. The representative function of the s^^ lattice cell is denoted as 9s'- 
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Summing over the monopole position according to Ref. ilTji], expanding the cosine function 
over the small fiuctuations in the fields x ^'^^ h, and integrating over these fields we get the 
monopole action: 
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where 



J^^l = Ss,s' - ml)b'^g, 
1 
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(11) 
(12) 



and Dma is the scalar propagator for a massive particle, {—df + m?')Dm{x) = S^^\x), with 
the Debye mass m = mo = A^^. Note that the lattice operators JF and Q are dimensionless 
quantities. 

In the infinite lattice case Eq. ()11|) can be rewritten as follows 
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where 



/i = h/X 
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The finite-volume expression for the monopole action can be obtained from Eq. (fT!?|) by the 
standard substitution: 
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(15) 



where L,- is the lattice size in i direction. 



In the infinite-volume case the lattice operator JF, ,,/ depends only on the dimensionless 
quantity /i, Eq. (fT^ . which is the ratio of the monopole size b and the Debye screening length, 
Eq.Q. The form of the operator JF is qualitatively different in the limits of small and large 
/i. Namely, the leading contribution to the monopole action is given by the mass (Coulomb) 
terms for small (large) lattice monopoles [2|: 

s,s' 

where -Dg ,,' is the inverse Laplacian on the lattice. Thus the Debye length A^ sets a scale for 
the lattice monopole size (or, better to say, for the size of the lattice cell) which characterizes 
different behavior of the monopole action. 

B. Squared monopole density in 3D 

The simplest quantity characterizing the lattice monopoles is the monopole density piattip) 
measured in the lattice units 

Plauih) = ^^{Y.\ks\). (17) 

where L is the lattice size in units of h. However, analytically it is more easier to calculate 
the density of the squared monopole charges, 

{k\h)) = ^^{Y.k\s))^ (18) 

which has a similar physical meaning to the monopole density. 

Using Eq.(j3I) the lattice density (J18p can be written in the continuum theory as follows: 

{k\h))= f d\ f dh{pix)p{y)), (19) 

where the lattice site s is fixed and the average is taken in the Coulomb gas of the magnetic 
monopoles described by the partition function Q. 

The correlator of the monopole densities, {p{x) p{y)), is well known from Ref. [l7|. In- 
troducing the source for the magnetic monopole density, J, Eq. ()19|) can be rewritten as 
follows: 



j=o 



Then we repeat the transformations in the previous Section which led us to Eqs. ()10|) . ()13|) . 
Integrating over quadratic fluctuations of the field x we get in the leading order 

(p(x) p{y)) = p [S^^'^x -y)- mlD^^ix - y)] . (21) 

Substituting Eq.(j2H) in Eq. lfTTHl and taking the integrals over the cell Cg we get 

{k') = pb' {J^-\,{p) , (22) 



where the inverse matrix JF, ,,/ is given by Eq. pi|) . Equation ()22j) estabhshes a direct relation 
between the density of the squared monopole charges and the monopole action ()lUp in a 
leading approximation of the dilute gas. 

The squared monopole charge satisfies the following: 



(^2^ ^,C,p\Db'-[l + {iXD/by)\ , ^ b:>\n, ^23) 



pb'-[l + C2p (b/XDY + O [ib/XDT)\ , 6 < Ad 

where Ci ~ 2.94 and C2 ~ 0.148 in the infinite lattice case. 

Equation ()23p can qualitatively be understood as follows. In the small b region the 
density of the squared lattice monopole charges is equal to the density of the continuum 
monopoles times the volume of the cell. This is natural, since the smaller volume of the 
lattice cell {Vol = b^) the smaller chance for two monopoles to be located at the same cell. 
Therefore each cell predominantly contains not more that one monopole, which leads to the 
relation k^ = \ks\ = 0,1. As a result we get (/c^) -^ Piatt{b) ^ pb^ in the limit 6 — > 0. 
In the large-6 region correlations between monopoles start to play a role. The monopoles 
separated from the boundary of the cell by a distance larger than Xd, do not contribute 
to (fc^). Consequently, the b^ proportionality for the random gas turns into Xob"^ in the 
Coulomb gas and we get {k^) ^ pXob'^- 

Finally, let us mention interesting relations between the density of the small- and large- 
sized monopoles and the coefficients in front of, respectively, the Coulomb terms and the 
mass terms of the monopole action, Eqs. (fTH |) .(P ^ : 

Cii>) = ^y b»Xn, and ^(6) = ^^ , 6 « A^ . (24) 

III. BLOCKING IN FOUR DIMENSIONS 

In four space-time dimensions the monopole trajectories are closed loops. Let us superim- 
pose a cubic lattice with the lattice spacing 6 on a particular configuration of the monopoles. 
Each of the (oriented) lattice 3D cells can be characterized by an integer magnetic charge 
it contains. Thus we can relate the continuum configuration of the monopoles to the lattice 
configuration. The three-dimensional cubes are defined as follows: 



Cs,^, = I b(^s^ - -) <x„ < b(^Su + -j for z/ 7^ /i ; and 



x^, = bSf,\ , (25) 



where s,^ is the dimensionless lattice coordinate of the lattice cube Cg^n and Xi, is the con- 
tinuum coordinate. The direction of the 3D cube in the 4D space is defined by the Lorentz 
index /i. 

The monopole charge Kc inside the lattice cube C^^^ is equal to the total charge of the 
continuum monopoles, k, which pass through this cube. Geometrically, the total monopole 
corresponds to the linking number between the cube C and the monopole trajectories, 
k (an illustration is presented in Fig. |2I). The mutual orientation of the cube and the 
monopole trajectory is obviously important. The corresponding mathematical expression 
for the monopole charge Kc inside the cube C is a generalization of the Gauss linking 




FIG. 2: Blocking of the continuum monopoles to the lattice in four dimensions. The lattice 
monopole charge is equal to the linking number of the monopole trajectory, k, with the surface of 
the three-dimensional cube C. 



number to the four dimensional space-time [5|: 



A-K^' 



d*x / d'^ye^,^ap^l^{x)K{y) 



\x-y\ 



(26) 



Here the function S^^(x) is the two-dimensional (5-function representing the boundary dC 
of the cube C. In general form it can be written as follows: 



^a(i\X) 



/d.,dr,^^«<')[.-i(f)l 
Jt. dTa On 



(27) 



where the four dimensional vector x{f) parameterizes the position of the two-dimensional 
surface S. The function V^'^' in Eq. ()2()|1 is the inverse Laplacian in four dimensions, 
d'pD'^^^x) = 6^^\x). It is obvious that the lattice currents Kg^^j^ are closed 



d'K = 



(28) 



due to the conservation of the continuum monopole charge, d^k^ = 0. In Eq. (j28j) the symbol 
d' denotes the backward derivative on the lattice. 

Following Ref. [5| we derive the lattice monopole action starting from a particular model 
for the monopole currents. We consider the model dual superconductor the partition function 
of which can be written as a sum over the monopole trajectories: 



Zmon = W^k j VB exp|- / d' 



X 



^A + ^k^{x)B^{x) 



^int[k) I 



(29) 



10 

where F^^ = d^B^ — dyB^ is the field stress tensor of the dual gauge field S^, and Sintik) is 
the action of the closed monopole currents fc, 

A;^(x) = ^dr^^^5(^)[a;-£(r)]. (30) 



Here the AD vector function x^(r) defines the trajectory of the monopole current. In Eq. ((2£ 
the integration is carried out over the dual gauge fields and over all possible monopole 
trajectories (the sum over disconnected parts of the monopole trajectories is also implicitly 
assumed) . 

The action in Eq. (|^^ contains three parts: the kinetic term for the dual gauge field, the 
interaction of the dual gauge field with the monopole current and the self-interaction of the 
monopole currents. The integration over the monopole trajectories gives the Lagrangian of 
the dual Abelian Higgs model J19J : 



Z^on cx Zdahm = JV^ JVB exp{-Jd'x [^F^^, + ^\{d^ + tB,)<^\' + V(<l>)} , (31) 

where $ is a complex monopole field. The self-interactions of the monopole trajectories 
described by the action Sint in Eq. (j^^ lead to the self-interaction of the monopole field 
$ described by the potential term V^($) in Eq. ()31|) . This model is nothing but the usual 
Ginzburg-Landau model written for the dual fields $ and B^. 

Similarly to the three dimensional case let us rewrite the dual superconductor model (|^ 
in terms of the lattice currents Kc, Eq. (J26|) . To this end we insert the unity, 



1= J2 UKKc-HdC,k)), (32) 

Kce^ c 

into the partition function ()29|) (here 5 represents the Kronecker symbol). Then we integrate 
the continuum degrees of freedom, k^ and B^, getting the partition function in terms of the 
lattice charges Kc- The simplest way to do so is to represent the product of the Kronecker 
symbols in Eq. ()32|) in terms of the integrals, 

oo 

E [n / d^c] exp{z Y^OcKc-ij d^x k^{x)B^{e- x)] , (33) 



Kc&Z C _oo C 



where 



B,ie;x) = ^ld%e,,apd,V^'\x-y) ^^cSf^(l/) • (34) 



To derive Eqs. (jH!^ . (jH^ from Eq. (j!?^ we have used relation (^Hjl . 
Substituting Eq. §^ into Eq. ^ we get: 



Z^on = ivk JVB J2 [n / d^c] exp{tJ2ecKc 

- j d\l^Fl, + ik,{x)[B,{x) + B,{e-x))]- SUk)} . (35) 
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One can see that the substitution of the unity ()33|) effectively shifts the gauge field in the 
interaction term with the nionopole current, B^ ^ B^ + B^. Therefore the integration over 
the nionopole trajectories, k^, in Eq. ()35|) is very similar to the integration which relates 
Eq. ^ and Eq. (jSH). Thus, we get: 



Zmon^ZoAHM = JV^ JVB ^ [j] / d^c] expjz J] 0^ i^C 
•^ •' Kc(^^ C J^ C 



- hTx 



42 M-^2 



d^ + i{b^{x) + B^{e- x))] $1' + V{^)\ } . (36) 



Next we rewrite the continuum dual superconductor model in terms of the lattice 
monopole currents, K: 



Zdahm — / . e 



^mon\J^ ) 



(37) 



where the monopole action is defined via the lattice Fourier transformation: 



"Ornon (-li ) 






(38) 



Here the action S{9) of the compact lattice fields 9 is expressed in terms of the dual Abelian 
Higgs model in the continuum: 



-m . 



P$ JVB exp{-/d 



X 



1 



V M-^2 



a,, + i(b^, + 5^(0))] $f + l^($)] } . (39) 



An exact integration over the monopole, $, and dual gauge gluon, _B^, fields in Eq. flH^ 
is impossible in a general case. Let us however consider the quadratic part of the monopole 
action. Neglecting the quantum fluctuations of the monopole field we work in a mean field 
approximation with respect to this field, $ — > ($): 



-s{0) . 



jVB exp{-/d^a;[— F^^, + Y (^/^ + ^/^(^))1 } 



Ag 



(40) 



where rj = \ ($) | is the monopole condensate. 

The Gaussian integration over the dual gauge field can be done explicitly. In momentum 
space the effective action (up to an irrelevant additive constant) reads as follows: 



^(«)4/(&^"'^-)%r"^ 



-p) 



(41) 



where B^{9,p) is related to the field Bfj_{9,x), given in Eq. (jBH), by a continuum Fourier 
transformation: 



^3 



P s,a 



(42) 
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with 



Q.(-) = n^^- (43) 



To get Eq. (jl^ from Eq. (jH^ we notice that 

^6^,„^ Sf^(x) = 9[^, Kf (x) , (44) 



where V^ is the characteristic function of the lattice cell Cg^fi- Namely, the characteristic 
function of the 3D cube with the lattice coordinate s^ and the direction a is 

V^iCs,a, x) = 5^,„ 6ix^ - bs^) n QiKh + 1/2) - x^) ■ Qix^ - bis^ - 1/2)) , (45) 

where B(x) is the Heaviside function. The Fourier transform of the function (j45|) is 

V^{C,,o.,p) = 5^,a b^ Qaipb) e-^'^P'^^ . (46) 

Integrating all variables but the lattice monopole field Kg,^ we get the quadratic monopole 
action: 

.'?..s' rv.ry' ' 



where 



•^.lU' = / (Cfi^^^pTT^^"^^^^"'^^^ '"'''"'' ' ^^^^ 



and 



^i = MBb. (49) 

In the fi —>■ oo limit the leading contribution to the operator J-" can be found explicitly: 

Sss',aa' = 21.2 -p ' ^aa'^s^^s'^ ' zZ ^Si^Sj^s^, , (50) 

cyclic 

where F = F(0, tt/y^l^^)? T^a\si_) is the three-dimensional Laplacian acting in a times- 
lice perpendicular to the direction d, 5s is the Kronecker symbol, Ag = P*^^^(s) is the 
one-dimensional Laplacian operator (double derivative), F(a,a;) is the incomplete gamma 
function and tjjy is an ultraviolet cutoff. 



IV. MONOPOLE DENSITY AND (MAGNETIC) DEBYE MASS IN 3D 

GLUODYNAMICS 

A. Technical details of numerical simulations 

In the next three sections we are discussing numerical results for the Abelian monopoles in 
the SU(2) gauge model. These monopoles obviously possess much more non-trivial dynamics 
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than the monopoles in the simplest case of the cQED. Nevertheless, we show below that in 
a certain limit the dynamics of the Abelian monopoles in the 3D SU(2) gauge model can 
be described by the Coulomb gas as in the CQED3 case. As for the 4D SU(2) model the 
Abelian monopoles in this case can be described by the dual superconductor model. 

We simulate numerically the pure SU(2) gauge model in three dimensions on 48^ lattice 
with the standard Wilson action 5* = — l/2X]p Tr f/p, where Up is the plaquette matrix 
constructed from the gauge link fields, f/;. To study the Abelian monopole dynamics we 
perform Abelian projection in the Maximally Abelian (MA) gauge jl5| which is defined by 
a maximization condition of the quantity R\U] = TTj2s,ii[Ufi{s)a3Ul^{s + /i)o"3]. 



max 

n 



/?[[/(^)] . (51) 



with respect to SU{2) gauge transformations, U -^ U^^^ = Q,'^Ui^. The gauge fixing condi- 
tion (j51|) is invariant under an Abelian subgroup of the SU{2) group. Thus the condition 
(jHH) corresponds to the partial gauge fixing, SU{2) — *• [/(I). 

After the MA gauge fixing, the Abelian {Ufj_{s)} and non-Abelian {U^{s)} link fields are 
separated U^{s) = C^(s)m^(s) where 



V/I - M')\' -£;M_] „-. /e'-.W 



^"W^i' .; ./rretrp)- ""W^lo e-io.i- («2) 



The vector fields C^(s) and Ufj,{s) transform like a charged matter and, respectively, a gauge 
field under the residual U(l) symmetry. Next we define a lattice monopole current (DeGrand- 
Toussaint monopole) [20|. Abelian plaquette variables 6^„{s) are written as 

V(^) = ^/^(^) +0,is + ix)- e^is + z>) - e,{s) , (-47r < e^,{s) < An). (53) 

It is decomposed into two terms using integer variables n^^^s): 

0^,y{s) = e^^{s) + 2^xn^^{s), (-tt < ^^^.(s) < tt). (54) 

Here 6^u{s) is interpreted as an electromagnetic flux through the plaquette and n^y{s) cor- 
responds to the number of Dirac string piercing the plaquette. The lattice monopole current 
is defined as 

k{s) = -tupadunpa{s + fi) . (55) 

In order to get the lattice density for the monopoles of various sizes, 6, we perform numer- 
ically the blockspin transformations for the lattice monopole charges. The original model is 
defined on the fine lattice with the lattice spacing a and after the blockspin transformation, 
the renormalized lattice spacing becomes h = na, where n is the number of steps of the 
blockspin transformations. The continuum limit is taken as the limit a ^ and n ^ 00 for 
a fixed physical scale b. 

The monopoles on the renormalized lattices ("extended monopoles", Ref. [7|) have the 
physical size b^. The charge of the n-blocked monopole is equal to the sum of the charges 
of the elementary lattice monopoles inside the n^ lattice cell: 

n-l 

fc(")(s) = J2 k^ns + ifi + ju + lp). 

i,j,l=0 
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For the sake of simplicity we omit below the superscript (n) while referring to the blocked 
currents. We perform the lattice blocking with the factors n = 1 ... 12. All dimensional 
quantities below are measured in units of the string tension a. The values of the string 



tension are taken from Ref. 2l|, |22 . 



In order to get rid of the ultraviolet artifacts we have removed the tightly-bound dipole 
pairs from all configurations using a simple numerical algorithm. Namely, we remove a 
magnetic dipole if it is made of a monopole and an anti-monopole which are touching each 
other [i.e., this means that the centers of the corresponding cubes are located at the distance 
smaller or equal than \/3a). Note, that we first apply this procedure to the elementary a^- 
monopoles, and only then we perform the blockspin transformations. Below we discuss the 
results obtained for the monopole ensembles with the artificial UV-dipoles removed. 



B. Parameters of the monopole gas 

In Figure Et^a) we show the density of the squared monopole charges (without the UV 
dipoles and normalized by the factor b"^) as a function of the scale b for various blocking 
factors, n. One can see that the 6-scaling violations are very small. As the blocking size 
b increases the slope of the ratio {k1)/b'^ decreases in a qualitative agreement with the 
prediction from the Coulomb gas model ((221) • 





. 


1 1 1 1 1 1 1 


, , , 


1 1 1 1 1 




■ No 


dipoles 




^^^^^ 


0.2 


- 






- 


^5 

A 
V 

0.1 


/ 






n=l 

n=2 
n=3 
n=4 
n=6 




/ 






n=8 
n=12 


n 


V 











b 



0.15 


¥ 






n=l 






^^ 


T 


n=2 




V 


•< 
A 


n=3 
n=4 


V 




X 


* 


n=6 






V 


■ 


n=8 


0.1 


- 


N 


• 


n=12 






No dipoles 


X. 


^ 


n nc 









^^^ 



b 



(a) 



(b) 



FIG. 3: The density of the squared monopole charges, {k'^), with the UV dipoles removed. The 
density is normalized (a) by 6^ and (b) by 6^. The fits by the function (|18() are shown by dashed 
lines for each value of the blocking step, n. 



According to the prediction coming from the Coulomb gas model fj23|) the ratio {kf) /b^ 
should tend to a constant as b becomes smaller. This behavior can indeed be seen from 
Figure Efb). Note that at small values of b the 6-scaling of the monopole density is violated. 
This scaling violation is not unexpected due to the presence of the lattice artifacts at the 
scale 6 ~ a. In order to get artifact-free results we will use below large-6 monopoles. 

The values of the parameters of the Coulomb gas model in the continuum limit, 
Eq. 0, can be obtained by fitting the numerical results for {kf) by the theoretical pre- 
diction (fTHj) . (fT3j) . Technically, for each value of the blocking step, n, we have a set of the 
data corresponding to different values of the lattice coupling /3, and, consequently, to differ- 
ent values oib = n ■ a[[3). Note that by fixing n we simultaneously fix the extension of the 
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coarse lattice, L/n, in units of h. The size of the coarse lattice enters in Eq. ()15p. We fit the 
set of the data for the fixed blocking step n. The best fit curves are shown in Figures Efa) 
and (b) by dashed lines. The quality of the fit is very good, x^ /d.o.f. ~ 1. 

The fits of the density provide us with the values of the continuum monopole density. 



p("\ and the Debye mass, M)^ , obtained for the fixed blocking n. These results are shown 
(in units of the string tension) in Figures Ufa) and (b), respectively. We expect to get the 
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FIG. 4: (a) The density of the continuum nionopoles, p, and (b) the Debye screening mass, M/), 
obtained with the help of the fits of the n~blocked squared monopole density by function ((TH|l. The 
large-n extrapolation (|56|) is shown by solid lines. 



artifact-free results in the limit of large 6, or, in our case, in the limit of large n. Thus, 
one may naturally expect that in the limit ra ^ oo the values of p^"-* and M^ converge 
to the physical values: lim„^oo O^'^^ = O^^, where O stands for either p or M^. We found 
that the dependence of both p and Md on the blocking size n can be approximated by the 
dependence 



Q{n) ^ Qph ^ ^Qj^g^ . 



n 



(56) 



at n > 2 according to Figures |ll^a,b). Using the extrapolation (j3Bj) we get the physical values 
for the monopole density p and the Debye screening mass Md coming from the Coulomb 
gas model (here and below we omit the superscript "ph" for the extrapolated values): 



p/a-"^' = 0.174(2) , Md /a^^' = 1.77(4) . 



1/2 



(57) 



The value of M^ may be treated as the "monopole contribution to the Debye screening 



mass 



J) 



In order to check whether the monopole dynamics can be described by the Coulomb gas 
model (j21) we construct the dimensionless quantity |3| 



C 



Mna 
P 



(58) 



which is known to be equal to eight {C'"'^ = 8) in the low density limit of the Coulomb gas 
model [l7| . In Figure \^ we plot our numerical result for C as a function of n. 
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FIG. 5: The same as in Figure |1] but for the ratio H58() . The dashed hue corresponds to the low 
density hmit of Coulomb Gas model 17], C = 8. 



Using the large-n extrapolation (jKH|l we get 



C = 10.1(11 



I.e. 



C 



1.26(3) 



(59) 



The quantity C is about 25% larger than that predicted by the Coulomb gas model in the 
low monopole density approximation, C£p = 8. The discrepancy is most likely explained by 
the invalidity of the assumption that the monopole density is low. Indeed, the low-density 
approach requires for the monopole density to be much lower than a natural scale for the 
density, g^ (remember that the coupling g has the dimensionality mass^^'^). The requirement 
p ^ g^ can equivalently be reformulated as p/M^ ^ 1, which means that the number of the 
monopoles in a unit Debye volume. Voir) = -^i) = ^d^i must be high. Taking the numerical 
values for p and Md from Eq. (fH7|) we get: plM^ ^ 0.03 ^ 1. Thus, the low-density 
assumption is not valid in the 3D SU(2) gluodynamics. However, the discrepancy of 25% 
observed in the quantity C, Eq. (j3^ . is a good signal that the Coulomb gas model may still 
provide us with the predictions valid up to the specified accuracy. 

One can compare our result for the monopole density, Eq. (j37jl . with the result obtained 



by Bornyakov and Grigorev in Ref. [23|, p^^ = 2"'^(1 ±0.02) g^. Using the result of Ref. J2lj . 
a/o" = 0.3353(18) g'^, we get the value p^'^ /a^/"^ = 0.207(5), which is close to our independent 
estimation in the continuum limit (J57j) : p/ p^'^ = 0.83(4). The result of Ref. J22| is about 
20% higher than our estimation for the monopole density. Thus, although the condition of 
the low monopole density approximation is strongly violated, the BFC method (based on 
the dilute gas approximation) gives the value of the monopole density which is consistent 
with other measurements. 

It is interesting to compare the result for the screening mass (j37j) with the lightest glueball 
mass measured in Refs. |2l|, |2^, Mo++ = 4.72(4) y/a. In the Abelian picture, the mass of 
the ground state glueball obtained with the help of the correlator. 



{F'MFL{R))= const, e 



-M. 



0++ 



R 



+ . 
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must be twice bigger than the Debye screening mass, 2Md/Mo++ = 1, where the Debye 
mass is given by the following correlator 

{F^uiO) F^,{R)) = const, e'^^^ + .... 

The comparison of our result (jKTj) with the result of Refs. J2l|,|22] gives 2Md/Mo++ = 0.75(4). 
The deviation is of the order of 25% similarly to case of the quantity C. 

Finally, let us compare our result for the monopole contribution to the Debye screening 
mass, in Eq. (jSTj), with the direct measurement of the Debye mass in 3D SU(2) gauge model 
made in Ref. 2J|, rriD j \fo = 1.39(9). The values agree with each other within the 25 



per cent: mct/mj^ = 1.27(11). Approximately the same accuracy is observed in the four- 
dimensional SU(2) gauge theory for the monopole contribution to the fundamental string 
tension hi 



C. Short summary 

The results of this Section indicate that the dynamics of the Abelian monopoles in the 
three-dimensional SU(2) gauge model can be described by the Coulomb gas model. Using a 
novel method called as the blocking of the monopoles from the continuum, we have calculated 
the monopole density and the Debye screening mass in the continuum using the numerical 
results for the (squared) monopole charge density. We have concluded that the Abelian 
monopole gas in the 3D SU(2) gluodynamics is not dilute. The self-consistency check of our 
results shows that the predictions of the Coulomb gas model for the monopole density and 
the Debye screening mass are consistent with the known data within the accuracy of 25%. 

V. STATIC MONOPOLES IN HIGH TEMPERATURE 4D GLUODYNAMICS 

A. Details of simulations 

Finite temperature system possesses a periodic boundary condition for time direction 
and the physical length in the time direction is limited to less than 1/T. In this case it is 
useful to introduce anisotropic lattices. In the space direction, we perform the blockspin 
transformation and the continuum limit is taken as a^ — > and Ug —>■ oo for a fixed physical 
scale b = n^Os. Here a^ is the lattice spacing in the space directions and Ug is the blockspin 
factor. In the time direction, the continuum limit is taken as at —>■ and Nt ^ oo for a fixed 
temperature T = l/{Ntat). Here aj is the lattice spacing in the time direction and A^^ is the 
number of lattice sites for the time direction. In general at ^ a^ (anisotropic lattice). After 
taking the continuum limit, we finally get the effective monopole action which depends on 
the physical scale b and the temperature T. 

The anisotropic Wilson action for pure four- dimensional SU{2) QCD is written as 

S = (3{- Y. P^J{s)+lY.P^d')}^ (60) 



■7 



s,i>j^4: s,i^4: 



P^,{s) = ^Tr[l - U^{s)U,{s + fi)Ul{s + v)Ul{s)] + h.c. (61) 
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The procedure to determine the relation between the lattice spacings a^, at and the param- 
eters /3, 7 is described in Ref. J26|. 

The monopole current is defined similarly to the three-dimensional case, 

k^is) = -e^upadunp„{s + fi) . (62) 

The monopole current satisfies the conservation law, d' k^{s) = 0. 

At a finite temperature the blockspin transformation of the spatial and temporal currents 
should be done separately [26]: 

71s — 1 nt — 1 

Kf,j^4{ss,S4) = Yl "Y k^^4^{nsSs + {ns-l)fi + ii> + jp,ntS4 + l), (63) 

ij=0 1=0 

Us-l 

Ki{ss,Si) = Yl h{nsSs + ifi + ji' + lp,ntS4 + {nt-l)), (64) 

i,j,l=0 

where Ug (n^) is the number of blocking steps in space (time) direction. 

We consider only the rit = 1 case since we are interested in high temperatures for which 
the monopoles are almost static. The lattice blocking is performed only in the spatial 
directions, Ug = 1 . . .8, and we study only the static components i^4 among the AD monopole 
currents i^^ (below we denote K4 as k.). At high temperature we disregard the spatial 
currents Ki since they are not interesting from the point of view of the long-range non- 
perturbative spatial physics. The size of the lattice monopoles is measured in terms of the 
zero temperature string tension, (Tt=o- 



B. Monopole action 

First, let us discuss the action for the static monopole currents. This action at high 
temperatures was found numerically in Ref. [26| using an inverse Monte-Carlo procedure. 
It turns out that the self-interaction of the temporal currents can be successfully described 
by the quadratic monopole action: 

Smon{k) =Yfi ^i(^) ' (^5) 

i 

where Si are two-point operators of the monopole currents corresponding to different sep- 
arations between the currents. The term Si corresponds to the zero distance between the 
monopoles, 5*2 corresponds to the unit distance etc. (see Ref. [26| for further details). The 
two-point coupling constants /, of the monopole action are shown in Figures. Efa,b) as a 
function of the distance between the lattice points. The numerical data corresponds to low- 
est, T = 1.6Tc, and highest available temperatures, T = 4.8Tc. The spatial spacings of the 
fine lattice, a^, range from a^ = .16cr~^/^ to a^ = .25cr~^/^. 

According to Eq. (J16j) the leading term in the monopole action for large lattice monopoles 
{b ^ Ad) must be proportional to the Coulomb interaction. 
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FIG. 6: Two-point coupling constants, fi, of the monopole action vs. the distance between the 
lattice points, r (in lattice units) for n^ = 6, various spatial spacings, a^, of the fine lattice. The 
temperature is (a) T = 1.6Tc and (b) T = 4.8Tc. The fits by the Coulomb interaction (|66j) are 
visualized by the dashed lines. 

To check this prediction we fit the coupling constants fi by the Coulomb interaction ()fifi|) 
treating Cc as a fitting parameter. The fits are visualized by the dashed lines in Fig- 
ures infa,b). As one can see from the figures, this one-parametric fit works almost perfectly. 
By fitting the action, we obtain the values of Cc for a range of the lattice monopole 
sizes, by/a = .96 . . . 1.5 and the temperatures, T = (1.6 . . . 4.8)Tc. According to Eq. lfTBj) the 
pre-Coulomb coefficient Cc{b,T) at sufficiently large monopole size b ^ Xn must depend 
on the lattice monopole size b as follows: 



Ccib,T) 



1 



R{T) 62 



where R is the product of the screening length and the monopole density 

R{T) = XniT) piT) . 



(67) 



(68) 



We present the data for the pre-Coulomb coefficient, Cc{b,T) and the corresponding 
one-parameter fits (|67p in Figure [7|^a). The fit is one-parametric with R being the fitting 
parameter. Again we observe that the agreement between the data for Cc and the fits is 
very good. We show the quantity R vs. temperature in Figure [7||^b). 



C. Monopole density 

An independent information about the monopole dynamics can be obtained from the 
behavior of the lattice monopole density at various lattice monopole sizes. According to 
Eq.(j2ni) the large-6 asymptotics of the quantity (/c^(6))/6^ can be used to extract the prod- 
uct of the screening length and the continuum monopole density R, Eq. ()68|) . We plot in 
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FIG. 7: (a) The pre-Coulomb coupling Cc and the fits of Cc by Eq. (|67|) for various temperatures, 
T. (b) The product of the screening length and the monopole density, Ea. (|HS)) . calculated from 
the monopole action (in units of the string tension) . 
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FIG. 8: The ratio {k'^{b))/b'^ vs. lattice monopole size, 6, at (a) T = 1.6 Tc and (b) T = 4.8 Tc, 



Figures |Hfa,b) the quantity {k^{b))/b'^ vs. the lattice monopole size b for lowest and highest 
available temperatures. 

Our theoretical expectations (j^Hj) indicate that the function {k'^[b))/b'^ must vanish at 
small monopole sizes and tend to constant at large b. This behavior can be observed in 
our numerical data, Figure |H1 The large-6 asymptotics of {k'^{b))/b'^ allows us to get the 
quantity ^] R in Eq. ^. 
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D. Check of Coulomb gas picture 

Let us denote Ract (Rp) the quantity R obtained from the behavior of the monopole action 
(density). From a numerical point of view these quantities are independent. Theoretically 
we expect that these quantities are equal. To check the self-consistency of our approach we 
plot the ratio of these quantities in FigurelHl^a). It is clearly seen that the ratio is independent 
of the temperature and very close (with 10%-15% deviations) to unity, as expected. 
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FIG. 9: (a) Check of self-consistency: the ratio of the quantities R, Ea. ()68() obtained from the 
lattice monopole action and density; (b) Check of the dilute Coulomb gas picture: quantities Csp, 
Eq. (J69|) . calculated from the action and density. 



A check of the validity of the Coulomb gas picture can be obtained with the help of the 
quantity 



Csp{T) 



a 



sp 



(T) 



a 



sp 



(T) 



Az,(T)p(T) RiT) 



(69) 



where asp is the spatial string tension. This quantity is similar to the one discussed in 
Eq. ()58j) in the previous Section. In the Abelian projection approach the spatial string 
tension should be saturated by the contributions from the static monopoles. In the dilute 
Coulomb gas of monopoles the string tension is [l7[: a = 8y/p/gM while the screening length 
is given by 0- These relations imply that in the dilute Coulomb gas of monopoles we should 
get Csp = 8. 

We use the results for the spatial string tension of Ref. [28[ in the high temperature SU{2) 
gluodynamics. It was found that for the temperatures higher than T ^ 2Tc the spatial string 
tension can be well described by the formula: aspiT) = 0.136(11) ■ gi^iT) -T^, where g^niT) 
is the four-dimensional SU{2) 2-loop running coupling constant, 



94DiT) 
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with the scale parameter A^- = 0.076(13) T^. Taking also into account the relation between 
the critical temperature and the zero-temperature string tension J29|, T^ = 0.69(2) ^/a we 
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calculate the quantity Cgp and plot it in Figure IHtb) as a function of the temperature, T. 
If the Coulomb picture works then Csp should be close to 8. From Figure El^b) we conclude 
that this is indeed the case at sufficiently high temperatures, T /Tc ^ 2.5. 

E. Short summary 

The main result of this Section is that temporal currents of the Abelian monopoles in 
the SU(2) gluodynamics at high temperatures can be described by the three dimensional 
Coulomb model with a good accuracy. This result indicates that the non-zero value of 
the three dimensional (spatial) string tension at high temperatures is due to the temporal 
Abelian monopoles. 



VI. MONOPOLE CONDENSATE IN 4D GLUODYNAMICS 

Finally, let us consider the SU(2) gluodynamics at zero temperature. The value of the 
monopole condensate 77 was previously estimated from the chromoelectric string analysis of 
Ref. |27(] to be 77 = 194(19) MeV. Below we determine the value of the monopole condensate 
from the effective monopole action. We skip a description of the numerical simulations since 
it is quite similar to the one discussed in previous Sections (we use the isotropic Wilson 
action for the gauge fields and fix the Maximal Abelian gauge). We mention only the 
explicit construction of the extended n^ monopoles 
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FIG. 10: The fits of the n = 6 monopole couplings by function (|5n)) . 

We get the quadratic monopole action using the inverse Monte-Carlo simulations. The 
definition of the couplings fi of the monopole action is quite similar to the three-dimensional 
case discussed in previous Sections. The couphngs are described in detail in Ref. |5|. We 
illustrate the success of the method showing the fitting of the couplings by the theoretical 
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prediction ()5Up in Figure ^| The best fit parameters obtained from tlie fits of different 
couplings /, are very close to each other. This fact provides a nice self-consistency test 
of our approach. The numerical value of the monopole condensate turns out to be n = 
243(42) MeV. This value is very close to the value r] = 194(19) MeV obtained in Ref. Hj 
using a completely different method. 

VII. CONCLUSIONS 

The BFC method together with numerical simulations turns out to be a useful tool to 
obtain non-perturbative information about the topological defects in the continuum limit. 
The application of this method to the Abelian monopoles in SU(2) gauge model gives rise 
to the following results: 

1. In the three dimensional SU(2) gluodynamics the Abelian monopoles can be described 
by the Coulomb gas model. The monopoles do not seem to be in the dilute gas regime. 
Nevertheless, the continuum values of the monopole density (p = 0.174(2) cr^/^) and 
the Debye screening mass {Md = 1.77(4) o"^/^) - obtained with the help of the dilute 
monopole gas model - are consistent within the accuracy of 25% with the known data 
obtained from independent measurements. 

2. In the four dimensional SU(2) gluodynamics the static Abelian monopoles can also 
be described by the Coulomb gas model at high enough temperatures, T ^ 2.5 Tc. 
The monopoles form the dilute gas. The spatial string tension - obtained in indepen- 
dent measurements - is consistent with the prediction of the monopole Coulomb gas 
model. In other words, in the continuum the spatial string tension is dominated by 
contributions from the static monopoles. 

3. In the four dimensional zero-temperature SU(2) gluodynamics the value of the 
monopole condensate, rj = 243(42) MeV, was obtained in the framework of the dual 
superconductor picture. This result is consistent with the result obtained previously 
by an independent analysis. 
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